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ABSTRACT 


The feasibility of utilizing a gaseous core nuclear 
reactor to provide high enthalpy, high pressure gas flow 
for simulating atmospheric re-entry conditions was inves- 
tigated. The test facility uses a mixture of nitrogen and 
uranium in a closed cycle with no attempt to contain the 
uranium fuel within the core. The primary purpose of the 
facility is to provide high enthalpy, high shear flows 
for testing re-entry materials and shapes. 


Investigated in this study were the effects of the 
nitrogen-uranium mixture on reactor criticality, nuclear 
contamination of the test model, protection of the reactor 
core and nozzle structure from imposed heat loads and operating 
limitations of the test facility. | 
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Title: Associate Professor of 
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SYMBOLS 
Avogadro's number, 6.023 x 107° 
molecules/gram mole 
area, eee eee 
atomic number, species 
buckling, Ben eneeeeie 
coefficient 
constant | 
diameter, centimeters 
modulus of elasticity, kilograms/millimeter, 
penne fone anche 
blowing eae 
transpiration cooling mass eco 
enthalpy, gram-caleries/gram 
earth gravitational constant, 4 x oe 
menere eeaonce 
neutron diffusion time, seconds 
diffusion length, centimeters 
length, centimeters 
molecular weight, grams/gram-mole 
number density, par enlelee cameras 
pressure, absolute atmospheres 
power, megawatts 
heat flux, gram-caleries/centimeter“second 
radius, centimeters | 
reflecter saving, centimeters 
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temperature, Kelvin 
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velocity, centimeters/second 


volume,centimeters 

work, megawatts 
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speed of sound, centimeter/second 
thickness,centimeters 

thermal absorptivity 

acceleration constant, a 7s0.7 centimeter) -eaorta 
coefficient of convective heat transfer, 

Bene ics cont ieea eae sec - ° Kelvin 

thermal conductivity, caleries/centimeter - Kelvin 
BOlezmann s Constant, ..3905e oes ergs/°Kelvin 
multiplication factor 

mass flow rate, grams/second 

pressure differential, atmospheres 
radius,centimeters 

time, seconds 

neueron velocity, meter/seconds 

axial distance from nozzle throat, centimeters 
altitude, nautical miles 

coefficient of linear expansion centimeter/ 
centimeter- Kelvin 4 

density proportionality cone eemae OL ALS, 

density exponential senttenta 1/7 .3 x 10°meters 
ratio of specific heats 

extrapolation distance, centimeters 

energy per fission, megawatt-sec/fission 
absorptivity 

mole fraction 


fission fragment fraction 
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Subscripts 
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angle, radians 

mean radiation absorption coefficients, cm 
decay constant 
mean free path, centimeters 

viscosity, gram/second-centimeter 

neutron fission multiplication factor 
density. Sans) cones 

stress, kilograms/millimeter, pound saeor Ce amare 
Stepfan-Boltzman Constant, so cums 1Ome erg/ 
are ~ °Kelvin - sec 

neutron flux, neut one cent ameter - second 
temperature-Mach number correction exponent 
angular velocity, radians/second 

macroscopic cross section, eenenmeee nein 
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half life,seconds 
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INTRODUCTION 


Many problems areas have been uncovered and investi- 
gated in our attempt to master space. One area is the 
re-entry of a ballistic, orbital or superorbital body 
into a planetery atmosphere and the delivery of an 


payload to the planet's surface. 


The kinetic and potential energy of the body must 
be dissipated to insure successful payload delivery. Atmos- 
pheric braking, with the transfer of energy from the high 
velocity re-entry body oo aa surrounding air, is often 


the optimum entry system. 


The body must be protected from the high heat loads 
due to friction and the rapid compression heating of the air 
immediately surrounding the body. Ablative materials are 
medeto both protect the body, by acting as am insulator, 
and to carry away the heat that is absorbed by the ejection 


Gmemass Of the ablative material. 


The ablative process, however, is very complex and 
attempts to analytically predict the behavior of new 
Materials and shapes have met with ieee success? Various 
attempts at simulating re-entry conditions have been made, 


but no existing facility can completely duplicate the 
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re-entry environment, although some aspects of the problem 


; aS 
can be simulated. 


Existing test facilities cannot duplicate the high 
enthalpy, high shear regime for a sufficiently long test 


time to be able to investigate this critical ablation region. 


This study will investigate the feasibility of a high 
enthalpy test facility which uses a fissionable gas reactors 
as a power source. The gas core reactor drives a closed 
cycle test system and provides the necessary energy for 


high enthalpy, high pressure gas flows of long duration. 


In establishing the feasibility of this system, the 
re-entry problem is first investigated to establish the 


power requirements of the system. 


Stagnation pressure and enthalpy are determined by 
the method outlined in Reference dee from these the 
power required to simulate an atmospheric entry in a ground 


test facility is determined. 


Then the reactor itself is studied. The gaseous core 
reactor utilized is similar in operating principle to 
the gaseous core nuclear rocket engine proposed as the 
next generation space ee However, as the 
proposed system incorporates a closed operating cycle, 
containment of the nuclear fuel does not present an obstacle 
as it does in the gaseous core nuclear rocket. Critical 
size, fuel concentration and fuel-nitrogen ratio, where 


nitrogen is the working fluid, are determined using the 





: : 9 
diffusion approximation and one group neutron theory. 
The contamination problems’ associated with the utilization 


of a nuclear medium are also investigated. 


The combined heat loads to the reactor and nozzle 
are examined. The heat loads in the nozzle are Ber omminedae 
and homogenous transpiration coon me to protect the nozzle 
1s examined. The use of a seeded carbon particle layer 
to protect the moderator from radiation heat loads is also 


; 2S 
examined. 


Finally the power requirements of auxiliary equipment 


to maintain steady state closed cycle operation are examined. 
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RE-ENTRY CONDITIONS 


{.l Pressure and Enthalpy 


To test the materials used in a re-entry vehicle it 
is necessary to know the velocity, the pressure and the 


heating involved in the re-entry. 


As a typical example consider a vehicle returning 
to earth with a superorbital velocity. An estimate of 
the velocity can be obtained by assuming it equal to the 
velocity needed to escape from a nonrotating earth, with 
no drag forces. The energy per unit mass to escape from 


the earth, acting as a central force field is 


12 & 
7 Ve R (1.1) 


and the velocity to escape is 


Ve = [2 2) 


when K is the earth gravitational Constant. 4x10" "Yes 


and R is the earth's radius,G.4 x10°m,. , the escape 





velocity £rom the Carth!’s surface ominous ole m/sec. 


This is the minimum for earth escape and is the mini- 
mum at which a superorbital body will approach the sensible 


atmosphere, commencing re-entry. 


Expressions for stagnation pressure and enthalpy, at 
the point of maximum pressure and also at: the point of 
maximum heating rate, may be found by using the velocity 
altitude history given in Reference] for an exponential 


atmosphere, 


2™ Sin Ge (1.3) 


-by 
-CrA Fo XE , 
BB 
V = Ve 


where Ve is the initial entry velocity, Cy is the coeffi- 
Stent Of drag of the body, m the mass of the body, O6¢e 
the entry angle and & and A are constants defined by | 


P = Pe a e7®% (1.4) 


The stagnation pressure on the body is approximated 
by Newtonian theory with the body angle considered small 
and the ambient pressure neglected. The stagnation pres- 


sure at any point in the entry is 


“Ps 
R ee 4 Vien exe (SA me (75) 


m Sin Ge 


, : AR 
For maximum stagnation pressure, set a4 =O and 


Sbptain 
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ee ces (pe os Pb CaS C,e°9-(p) =O (1.6) 


ay 
where 
Ce = P. « Ve 
CoS ioe: CoA & Do 
m 2 Be Bb 


C, and Cz are fixed by the atmospheric model 


Ballistic parameter, #25. and initial entry velocity; 


and remain constant during re-entry. 


The altitude for maximum stagnation pressure during 


re-entry is given by 


lo ee — — (ie) 


Since the velocity for maximum stagnation pressure is 


given by 
-3 
V=Ne€& tt -@) 


and associated density is 


e = msineey (1.9) 
Sof 


the maximum stagnation pressure is 


Dp a msin Be BNle (1520) 
MAX Co A e 
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and the stagnation enthalpy is 
2 
h + % (cll ita) 


Or since the static enthaply is very much less than dynamic 


enthalpy, 


2 
ya Ve (1.115) 


Similarly, expressions are obtained for stagnation pressure 
and enthalpy at the point of maximum rate of convection 
heating of the stagnation point, where stagnation point 
heating rate is given by* 


5 (11D 
= V/ 


where Ce is the stagnation point coefficient of heat trans- 
fr, a function of body configuration, G~ is the radius of 
curvature at the stagnation point, and © and V are functions 
Sere altitude. To obtain the point of maximum stagnation point 
heating rate, the previously derived expressions for SP 


and V are substituted in the above. 


d / dt. 
For the maximum set —( 977]/?Oobtaining 


dy 
aH. = MSinGe Ver 
cit Imam “4 coe ee (eS) 


which occurs at an altitude of 


x in GR | (1.14) 


™M BinBe p 








and velocity of 


V=Ve e es) 


The associated stagnation pressure is 


Mm Sin O Van | 
ee ee si (1.16) 


CoA Dae 
and stagnation enthalpy is given by 


7 Var (1.17) 
2e% 


Figure } shows the extremes of stagnation pressure 
verses stagnation enthalpy which are encountered during 
an atmospheric re-entry. The ballistic parameter, Somes f 
1s a free parameter in determining the envelope. ‘Two 
envelopes are shown, one for Ve=li2xlO m/sec , corresponding 
to velocity entering the sensible atmosphere equal to 
escape velocity. The other is for Ve=0.e8* 10+ m /sec, 
corresponding roughly to !00n."™. orbital velocity. It 
can be readily seen in Figure i that: 
1) stagnation enthalpy at maximum heating rate and 
Maximum stagnation pressure is independant of 
the ballistic parameter, but increases with in- 


creasing initial entry velocity. 


ii)stagnation pressure increases directly as entry 


velocity and directly with ballistic parameter. 


i.2 Power Requirements for Simulation 


The stagnation pressure on the test model is fixed 





by re-entry conditions but the chamber pressure is limited 
by structural considerations. The maximum chamber pressure 
is set at 1000 atmospheres for the analysis. It is desir- 
able to test at the lowest possible Mach numbers consistent 
with representing hypersonic pressure and heating rate dis- 
tributions. For this reason, and to minimize nozzle losses 
and cooling requirements, a test Mach number, Me , Of -3..0 


is used in this analysis. 


For isentropic frozen flow through the nozzle, the 


power required per unit throat area is 


RH. } 
SD Gen / YANG) 2 





If Equations 1.16 and 1.17 are used for cc and 
meee cespectively, the power required per unit throat area 
to simulate the conditions at the point of maximum heating 


mace 1s 


(1.19) 





These power requirements are plotted in Figure 2 against 
initial entry velocity with the ballistic parameter as a 
free parameter. Reference 14 was used to obtain the charac-— 


teristics of nitrogen,which was taken as the test gas. 


Limiting the chamber pressure to 1000 atmospheres 
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established one limit on the re-entry profiles that can 
be simulated. Another is-set by limiting the chamber 
temperature to 10,000°k, this upper limit is set after 
a consideration of the heat transfer loads experienced 


in the various test facility components. 


ies Test Time 


The main purpose of this facility is to examine the 
behavior of re-entry materials and shapes under the high 
heat and shear loads of an atmospheric re-entry.. Such 
materials are ablative in nature, that is, they absorb 
heat by increasing in temperature and changing in chemical 
Or physical state. The heat is carried away from the 
surface by a loss of mass. The departing mass also blocks 
part of the convective heat transfer to the remaining mate- 
pre The tunnel running times are set by the time required 
for the rate of ablation to reach equilibrium. This time 
is a function of the ablative acto and for typical 
ablative materials , such as, graphite and glass reinforced 
phenolics, the ablative equilibium times is extended 


to 200 seconde 


1.4 Existing Facilities 


Shock tubes are used to generate higher temperatures 
and pressures to simulate re-entry conditions. The Cornell 
Aero ae foot hypersonic tunnel can generate 
stagnation temperatures and pressures up to 4200°K and 
2000 atmospheres respectively. The test time for these 


pressures however is only 4 micro-seconds. This tunnel can 
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be upgraded to 10,000°K and 2000 atmospheres but still has 
a run time far too short to produce the time history of 
high temperatures and pressures at which the ablative 


materials fail. 


The test time for a shock tube may be increased by 
using several tubes fired in a Gatling gun fashion. This 
however only increases ae time to seconds, developing 
a jet power of 25, ooo HP” , about four orders of magnitude 


too small for re-entry requirements. 


High enthalpy and pressure flow is also possible 
ferbizing a hot shot tunnel powered by an electrical 
discharge arc. This type of tunnel requires large amounts 
of stored electrical energy and is presently limited to 


a few hundredths of a second Sneretonea 


Although these facilities can approach the desired 
levels of enthalpy and pressure, they cannot maintain these 
conditions for a sufficient lenght of time. They are severely 
power limited. As an example, the power required, from 
Figure 2, to simulate the re-entry of a vehicle whose ball- 
istic parameter is 2000 grams/ em* and initial entry 
velocity is 8000 meters/second, is 30 megawatts per ae of 
throat area. This is about ten percent of the total output 
of the Grand Coulee Dam for each square centimeter of nozzle 


throat area. 


Plasma arc tunnels provide high enthalpy but low 
pressure flows and are currently being used to test re-entry 
materials. These facilities are also power limited and operate 


at low stagnation pressures and use argon as a working fluid 
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to achieve long duration tests with comparatively low 


power inputs. 


Because of these high power requirements, a gaseous 
core reactor is suggested to drive the test facility, the 
gaseous core because its operating temperature is limited 
only by heat loads to the containing structure and can 


operate at the desired temperature, 10,000°K. 
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CHAPTER FE 


TE SimeraAC Tiel 


2.1 Gaseous Core Reactor Requirement 


The power required to produce long duration high 
pressure ,jhigh enthalpy gas flows appears to be technically 
achievable only with nuclear power. Conventional combus- 
tion is limited by the recombination of combustion products, 
which absorb energy and result in an equlibrium tempera- 
ture far below the required stagnation temperature. [In 
atomic reactors, solid and liquid core operating tempera- 
tures have upper limits, the former is the melting poOmnt 
of the nuclear core and, or, structure and the later is 
boiling point of the ligquidized fuel, and, or,melting point 


@eetehe structure. 


A gaseous core nuclear reactor is capable of steady 
state operation at very high temperatures and pressures. 
The operating temperature of a nuclear reactor utilizing 
a fissionable gas core is limited only by the heat load 
eiae the hot gas imposes on the container. These heat loads 
can be controlled to some degree by a combination of con- 
vective and transpiration cooling and by injection of a 
thermal radiation absorber to block radiation heat transfer 
to the wall. The gaseous core reactor in this study is 


capable of providing the energy necessary to run a test 
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at the desired stagnation pressures and enthalpy, that is, 


7,100 megawatts. 


mee Gas Cycle 


To simulate the effects of atmospheric heating, nitrogen 
is used as the working medium. Uranium 235, initially a | 
fine power, is the nuclear fuel. Biological hazards 
and economic necessity make it mandatory that the uranium 
fuel be contained within the aoe of a gas core nuclear 
Pee. As the proposed system is a closed system, these 
considerations do not apply. This is not to Say, a priors, 
that fuel containment in the core properis not required. 
Presence of uranium in the nitrogen flow complicates the 
analysis and interpretation of experimental data. If the 
ratio of mass flow of uranium to nitrogen is limited to 
a maximum value of ten percent, or about one percent mole 
ratio, it should permit the assumption that the flow in 
the nozzle and test section behave as a pure nitrogen flow. 
This limit is imposed as the.maximum allowable uranium 
content for thermodynamically simulating the effects of 


atmospheric entry. 


The presence of fissionable gas in the closed cycle 
causes induced radioactivity of the system components. 
State of the art shielding and radioactive materials handling 
peennicques are assumed sufficient to provide necessary 
biological protection. As this facility would be very 
expensive to construct,it is necessary that the induced 
radioactivity of the model be within biological limits 


for examination of the model within a short time after 








the test and associated radioactive exposure. 


The radioactive damage or effects to the model must 
be slight enough so that these associated effects do not 
mask or override the aerodynamic effects which are being 


“investigated. 


These nuclear effects on the model are investigated 
to determine their influence on the design of the over-all 
system, and to determine if containment, greater than that 
implied by the thermodynamic limit of ten percent .uraniun, 


1s required.’ 


Chamber pressure and temperature are set at 1000 
atmospheres and 10,000°K. A portion of the heat load to 
the moderator is absorbed by the flow of nitrogen entering 
the core, the remainder is removed by an external water flow 
enroute to another system component, the water diffusor. A 
seeded layer of carbon is used to provide a radiation barrier 


which protects the moderator from radiation heat loads. 


The thermal energy of the hot gas in the reactor core | 
is converted to kinetic energy in a convergent-divergent 
nozzle. The nozzle is cooled by conventional backside 
cooling and transpiration cooling, both utilizing a secondary 


flow of high pressure, but low enthalpy nitrogen. 


A schematic of the test facility is given in Figure 3. 
A temperature - enthalpy diagram of the nitrogen cycle is 


given in Figure 4. 
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After the test section, the secondary nitrogen flow 
is vented into the main flow, reducing the temperature 
slightly. The stagnation temperatures are still in the 
wicinity of 10,000°K and prevent the use of conventional 


turbo-expansion devices or heat exchangers at this point. 


A water spray diffusor is utilized to slow the flow 
and lower the stagnation pressure and temperature of flow 
neo 1000°K and 100 atmospheres. The water spray is sufficient 
tO cool the walls of the diffusor. 


A centrifugal separator is used to seperate Ehe 
condensed uranium particles from the flow of superheated 


steam and nitrogen leaving the water diffusor. 


At this point the flow is basically free of uranium 
and at a temperature and pressure at which conventional 
turbo-machinery can be employed. The flow has 2,500 megawatts 
of available power, if expanded only to the saturation point 
of the water. After this energy is removed from the stream, 
the water is condensed and separated. The nitrogen, at 
0.24 atmospheres and 370 K,is cooled further in a conven- 
tional heat exchanger and pumped back to 1000 atmospheres 
using two multistage centrifugal flow compressors with 
interpump cooling to maintain temperatures below 1000°K. 
The power required for the pump work, 880 megawatts, is 


provided by the aforementioned turbine. 


The pressure of the water is raised to 120 atmospheres 
in another multistage centrifugal pump, requiring 10 mega- 


watts of the above turbine work. 
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The temperature - enthalpy diagram of the water cycle 


is given in Figure 5. 


2.3 Component Interface Problems 


To be successful as a system, the components of the 
Facility must be mutually compatible. In the following 
sections, the operating characteristics of the reactor 
and the nozzle, which are the main components of the 
system, are examined. The reactor is examined to determine 
the effects of various operating pressures and fuel-nitrogen 


concentration on the critical size and fuel concentration. 


Contamination of the model is examined to determine 
mente imposes a constraint on the fuel—-nitrogen mixture 


am the core. 


Heat transfer in the reactor and nozzle are also 
examined to determine limits on operating temperature, 


power and component size. 
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CHAPTER 11a 


GASEOUS CORE NEUTRONICS 


3.1 General 


Reference 6 shows that the gaseous core nuclear 
reactor is theoretically feasible. It shows that one 
can obtain a reasonable estimate of reactor critical 
maaius and critical fuel concentration by considering 
the reactor core to be a homogeneous mixture of fuel 
and moderator or a reflected sphere of fuel alone. The 
effects of containing the fuel with a gaseous annulus 
of hydrogen are also investigated as were the effects of 
multiple cavities. Reference 6 quantitatively shows that 
a gaseous core reactor can be designed with either simple 


Or multiple cavity geometry. 


In this study, the reactor will be considered to be 
a homogeneous mixture of pure uranium 235 and Ne with 
a bare spherical geometry. The effects of reflecting 
the above bare sphere are also considered. One group 
neutron theory is utilized, that is, the neutrons are 
considered to be in thermal equilibrium with the reflector- 


moderator. 


No attempt is made to separate the fuel from the 


MWerking fluid, nitrogen. Criticality calculations 
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are made with the ratio of number density of nitrogen 

atoms to number density of fuel atoms, Nu /No / 

aS a parameter. This gives one measure of contamination 
of the flow in the model test section. From the fluid 
dynamic point of view, the flow is treated as behaving 

as pure nitrogen, and the number density ratio gives a 
measure Of how valid this approximation is. The model 

ts also subjected to radioactive contamination, the 

number density ratio gives a partial measure of this effect; 
partial in that the uranium itself in the flow is not 


the only source of model contamination. 


Smee Criticality Condition : 


As a first approximation, the reactor can be considered 
to be a bare homogenous mixture of pure uranium 235 and 
nitrogen. Criticality is determined by considering all 
of the neutrons in the core to be thermal, that is, there 
is no loss of neutrons in the process of slowing down 
from the high energy state of production in the fission 
process, to capture at a lower thermal energy state. For 
steady state operation, the conservation of neutrons requires 
Bee the rate of production of neutrons in a control volume 
be equal to the rate of loss of neutrons due to absorption 


and leakage from the volume. 


The rate of production of neutrons per unit volume 


is given by 


production = FY = Sal 
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where v is the average number of neieeone 
produced per thermal fission 
DB is the neutron flux 
and Zig 45 the macroscopic £issionserose@-caumron 


of the core material. 


The neutron leakage is given by the diffusion approxi- 


Mation as 


> : 
Leakage = -Dg Zp | (322) 

where y* is the La Placian operator 

and D is the diffusion coefficient, 


1/3 of the transport mean free path 


in the “reactor 


The absorption term is similar to the fission tern, 


Absorption = Bla +&.) (5) 


here S is the macroscopic absorption cross section 


Wrieeie sore. 


Equations 3.1,3.2 and 3.3 are combined and 
the conservation of neutrons is expressed in differential 


form by 


Be oS + (Voi. 2a = © (3.4) 


Or if the material buckling, Beis deeineae te 


B32 = Y(Z;-fa)/ D 3.5) 
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the conservation of neutrons becomes 


Vv" +B gd =0 (Gao) 


. For a spherical reactor, with the boundary conditions 
that the flux goes to zero at the radius of the core and 
the gradient of the flux is zero at the center of the 
sphere, Equation 3.6 can be solved for the geometric 


buckling of the reactor, 


By = = ° (3.7) 
The criticality condition is the matching of the 
geometric buckling which satisfies Equation 3.6 with 
the material buckling of Equation 3.5. . The diffusion 


length is defined as 
rc iD 
ee = iS we (358) 
and the infinite multiplication factor as 
Re = » i | (3.9) 
Zia 
With this notation, the criticality requirement is 


(We) = Rest (3.10) 


The criticality requirement for the system is that 
the geometric buckling be equal to the material buckling. 
This means that the rate of neutron creation due to fissions 
is equal to the rate of destruction of neutrons due to | 
absorption, that is, fission and nonfission capture, and 


leakage from the critical region. 


For a nitrogen and pure uranium 235 system, the number 


of neutrons released per thermal fission of uranium 235 is 1s 
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v= 2.4G fevcons 
wa “235 


The interaction cross sections per atom for uranium, 


sub-subscript wu, and nitrogen, sub-subscript w ,» are S 


oy 580 barns | Gs. = 10 barns 


680 barns Gs,, = 10 barns 


by 


es: 


“ 


Sa 
Ga, =1-88 barns 


where the subscripts § denotes fission, 
a. denotes absorption, 


s denotes scatter. 


With these cross sections, the criticality equation 
can be used to determine the number density ratio of 
nitrogen atoms to uranium atoms required for criticality 
as a function of the core radius and the number density 


Of uranium 235 atoms, 


Nu oy 199 a J AX\O' — a (3.1) 


Using this relation the critical radius is found as 





a function of uranium 235 concentration with Na /Nu 
aS a parameter. The results are given in Figure 6 with 
the mass flow ratio of nitrogen to uranium shown in lieu 


of the number density ratio. 


Examination of the criticality equation above shows 
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that as Nua Varies decreases below 199, RN must increase. 
Therefore, Nn / Nu = 199 gives the minimum reactor size 
for a given concentration, and the associated mass flow 
ratio is 15.4. This is due to the scattering contribu- 
mons of the nitrogen atoms in the core reducing the 
transport mean free path of the neutrons. This is most 


readily seen if the criticality condition is expressed by 


QR = TY * 3(I- 4018) | . 
NJ ( \Y Th - Vauw -o Gc. )f Osu + Nahe Gsu (3.12) 


where the macroscopic cross sections are replaced by 
the number density times the microscopic cross section. 
The denominator has a maximum value when Na / No is 
199. For higher or lower number density ratios, the 


@reieical radius increases. 


As used here RK is not the actual radius of the 
homogenous sphere, but rather that radius plus the extrapolation 
distance. The extrapolation distance is the distance 
at which the neutron flux, which has a finite value at 
the physical boundary of the sphere, goes to Zero. Reference 


9 gives the extrapolation distance as 


3=2\3D (amie) 


The physical radius of the bare sphere is 


r= R- 213D ae 
This critical size and uranium concentration so 
obtained is only an approximation and is valid only-in 
determining the effect of nitrogen poisoning of the reactor 
On critical size and critical uranium concentration. The 


actual reactor will be operating at elevated temperatures 
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which will effect the interaction cross section,which 
generally vary inversely with neutron energy. A higher 
equilibruim temperature means higher neutron energies 
and, in general, lower interaction cross sections. Also 
the actual reactor will not be bare but must be enclosed 


in some way. 


The reactor core is in fact contained in a pressure 
vessel, for structural reasons the walls of which must be 
maintained at a low temperature relative to the reactor 
core. Lining the pressure vessel with beryllium, a 
moderator, which must also be cooled for structural reasons, 
provides a means of decreasing feats from the core and 


a means of providing low energy neutrons to the core. 


Nitrogen, being heavier than berylliun, some not 
absorb as much of the neutrons energy in a collision, 
therefore more nitrogen collisions are necessary to thermalize 
a neutron. Also nitrogen, even at very high pressures, is 
still: very much less dense than the solid beryllium moderator 
and at the higher core temperatures, the scattering cross 
section for nitrogen is less than that for the cold beryllium 
moderator. It is reasonable £o assume then, that the 
neutrons are inthermal equilibruim with the cool reflector 


and not with the hotter nitrogen-uranium mixture. 


The effects of the reflector-moderator on core radius 
can be found from a one velocity, two region consideration. 
As before in the core, the conservation of neutrons is 


given by Equation 3.6. 








ae 


In the reflector, there is) nowtisstemasleoematerial 


and the conservation of neutrons gives 
é = 
D- \ oD. a Tip PD. a O (32515) 
where subscripts f and ¢ stand for reflector and core 


respectively. 


With the inverse of the square of the reflector 


diffusion length, L,, defined as 


Ko = & 


the neutron conservation equations in the core and 


en | (3.16) 


reflector is obtained by dividing Equation 3.15 by Dr, 


TH +B @ =0 (Grol 


and 


we D, oan hae De = (3.18) 


= 


This set of differential equations can be solved 


for a criticality condition with the boundary conditions 
O< Be(r = 0) < o (3.19a) 


Belr=Re) = Prlr=Re) (3-19b) 


DeVPe| = D-VGZ-| og (3.19c) 


and 


Dr(Retde +5) =e (3.194) 
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where rae is the reflector thickness 


and R is the critical radius Of she eee mleetecd core. 


The solution of Equation 3.17 and Saks with 
boundary conditions of 3.19 gives an expression for 
the criticality ofa reflected spherical reactor in the 


transcendental form, 


D (ke Re coth ke(de +3 )+ 1) = DelI-R-Be coh eRe) (3-20) 


The reflected core radius can be expressed in terms 


Of the bare core radius, & 2 and the reflector saving,9%, 
as 
ei 7tt - 
R : (3.21) 
This value is substituted into Equation 3.20 and 


a transcendental expression for reflector saving is 


obtained, 
col 7S — ce © ke cae. sire (\- Bs ihe ome (3.22) 
T me --S) 


Equation 3.22 is solved graphically for the reflector 
Saving at each known r, then the reflected core radius, R,, is 


obtained by Equation 3.21 


The resulting critical radii are plotted in Figure .7 
against critical uranium concentrations for a 100cm 
beryllium reflector. Shown is the reflected core radius, 
with the mass flow ratio of nitrogen to uranium as a free 


Parameter. Also shown are lines of constant pressure for 
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an operating temperature of 10,000°K. 


Heat transfer in the reactor core, coupled with 
practical limits on the capacities of the system's pumps 


and turbines, lead to the choice of a minimum core radius. 


Structural: consideration dictates a low operating 
pressure, but this is not compatible with the desire 
to simulate the high stagnation pressures of a re-entry. 
Also, the problems associated with heat transfer in the 
core require high mass flows and low temperatures. For 
these reasons the operating pressure is set at the maximum 
possible with projected structural technology, that is, 


1000 atmospheres. 


While it is also desired that the mass flow rate of 
nitrogen to uranium be maximized so that the test facility 
Simulate atmospheric conditions, the requirements for a 
minimum core size dominate in this case. It is also noted 
that the critical core radius increases rapidly with a 
slight increase in mass flow ratio from the minimum of 15.4 

For these reasons the core radius is set at the minimum 
possible with a pressure limit of 1000 atmospheres. This 
radius is 170 centimeters, occuring at the mass flow ratio 


of 15.4 
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CHAPTER IV 


MODEL NUCLEAR CONTAMINATION 


4.1 Radiation Damage 


In the gaseous core nuclear reactor used as a power 
source for this system, there is no provision made in 
the core itself for the separation of fuel and fission 
fragments from the working fluid. This makes the design 
of the reactor much more simple then if some hydrodynamic 
means of separation of fuel from the working fluid is 
utilized. As the system is closed there is no problem 
of venting radioactive material to the atmosphere ,relieving 
one of the major constraints that faces the use of a gas 


core nuclear rocket. 


The presence of radioactive material in the fluid 
stream causes contamination of the model being tested. 
This contamination is examined to determine if the structural 
damage to the model will sSignicantly influence or mask 
the aerodynamic and thermal effects being tested. As the 
model must be examined after the test, the extent of radio- 


active contamination as a biological hazard also is examined. 


Radiation damage is caused by neutron bombardment 
and interaction, alphaand beta bombardment and gamma 


madlation. 
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The radiations may be divided into two groups. The 


light group consisting of betas, gammas and other electro- 
magnetic radiations, and the heavy consisting of neutrons, 
fission fragments, accelerated ions, and alphas. As a 
general rule the damage from the heavy particles is more 
severe that that from the light panbveileemaane The heavy 
particles displace the atoms of a solid crystal lattice 
from their normal position. In general this will result 
in an increased ultimate strength, a decrease in elongation 
and a reduction in area of the material, also eee 
hardness of the material increases and there is a decrease 
in impact strength. These effects are minimized if the 
neutron bombardment occurs at elevated temperature, as 


Hemthne case in this test facility. - 


Graphite, a prime candidate for a re-entry material, 
incurs crystallite damage during irradiation. Hardness 
and strength are increased while thermal and electrical 
conductivities decrease. Graphite shapes change in gross 
dimension during radiation. -All of these effects decrease 
in magnitude when the temperature of the graphite during 
radiation is increased, with almost no effect when graphite 
temperatures are Seve 1000°c. Also the integrated flux 
levels at which these effects begin to be noticable is 
high, of the order of tO newerens per square centimeter. ah 
As this test facility will provide very high model tempera- 
tures and neutron fluxes at the model of the order of ice 


neutrons/ cm* sec, it is expected that irradiation damage 


will have negligible effect on the test results. 
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4.2 Biological Hazards 


In determining the biological radiation hazard due 


to model contamination the following assumptions are 


made: 


ee 


Phe 


radioactive particles, such as, fuel, fission 
fragments do not adhere to the model upon 


impact, 


contamination is due to neutron-model 
interaction causing a transmutation to 


a possibly unstable isotope, 


1iil.contamination is acceptable if the model 


can be handled with minimum protection within 


a short time of the test run, 


iv. acceptable radiation levels from the model 


are 7.5 x ay rem/hour as recommended by the 


US National Committee on Radiation Protection. 


4.2aContamination Environment 


Neutron density at the model has three sources; 


direct streaming from the reactor core, neutron flux 


from the fissionable mass in the test section, and secondary 


emission from fission fragments. 


Direct streaming of neutrons from the core can be 


measured by noting the attenuation of the neutron flux 


by the gases in the nozzle between the model and the core. 








pl 


The differential change in neutron flux due to absorption 


by the gas is 
eg =a (4.1) 


The macroscopic absorption cross section of the gas 
varies with position in the nozzle as the density of 
the gas decreases. Assuming this variation is linear 
with axial nozzle position, X ,'and that the density 
goes to zero at the model, Equation 4.1 becomes 


Ax 


so = SZ. (1-45) dx oe) 


When integrated and evaluated at the model, the 


flux at the model is given by 
-YY22ra> 


= e 
J. - D. © (4.3) 
where L. is the axial position of the model and Aq 


is the absorption mean free path 


The axial position of the model, lu ysl meee 
by heat transfer considerations in the nozzle and inlet 
and must be minimized. Then lL is the order of the 
required nozzle length to expand the flow to Mach 3, 


approximately 30 cm, for the throat size set in Section 6.1. 


. For the chamber condition of this study the mean 
free path for absorption in the chamber is about five 
times this value, then the free streaming flux from the 


core to the model is approximated by 


Ds =O9 De (4.4) 
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The neutron flux from the mass of fissionable material 
in the test section is found by considering the change in 
flux as the flow goes from the critical region of the 


eore to the subcritical test section. 


The time of flight from the core to the test section 
is of the order of the length divided by the speed of 


sound at the nozzle throat, 
ty — Li /c* (4.5) 


If the test section is considered to behave as 
an infinite cylinder with constant cross section, the 
Critical buckling is a function of the geometry and is 


related to the test section radius by 


2 a 280s )* (4.6) 


The effective multiplication factor in the test region 


is ‘ ; 
ee es 
Kes = 2.405 D 
+ (2405 Re ad ) Cove a: 


where \K,, is defined by Equation 3.9 


With the assumption that the number of neutron generations 
is given by the time of fletqine yan divided by the average 
neutron Lifetime, Aa /Un , the neutron flux at the model, due 
to the reactivity of the fissionable gas in the test section, 
is given by 


_ ts (Regs-n (4.8) 
Din = Do Er 


Ro. Ret 
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Equations A.5S5 Ee haouG la 4.3 are evaluated for 


the conditions in this facility, 


_~ 4 
t = &x/O seconds 


keg = -9002 
Hep- to 
‘heey SOOm 


to obtain the flux at the model due to the reactivity of 


the test section, 50 


—| 
Ps - DHE - (4.9) 


As Psc is very much smaller than g¢s : ite 


is considered negligible in this analysis. 


Secondary emission of neutrons by fission fragments 
in the vicinity of the model will also GCONELICULeoeLoneme 
neutrons bombardment of the model. Given a fission fragment, 
mee «with a half life, ty a , and a yield fraction 
per fission of V[, ,from the definition of half life, the 





radioactive decay constant , XD» Pees 
2 
Tyo 0 


Mie rate of creation of atoms of fragment A is 


eo = Va Za Po (4.11) 


where 3 Gis the number of fissions per unit volume per 
unit time. 


The rate of decay is 
any 2 2a (4.12) 
dt 


mags in the reactor, fission products, A Toc ap at 


the rate 


\ 
cf = Ye 2a, oe — AKA (4.13) 
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and concentration of species A at time ewes 


AQ = x (me ZA - {Yn a, Do - -asAete* ad (4.14) 


With the initial condition that the initial concentration 

of species A is zero, the concentration of fragments leaving 
the core is given by Alto), where te is the dwell time 
of fluid in the core, | | 

poe Se Vere (4.15) 
> vr) 

At tp =lsecond, the concentration of A leaving the 


reactor is given by 


A()) = aioe (1-e™ (4.16) 


The rate of decay in the nozzle is given by Equation 4.12 
evaluating this at the test section with a time of flight, 
er, from the core, the rate of decay of fragment A 


at the model is 


dA _ Na Bi Ho a De Aw) rate (4.17) 
ct - ag me 

If it is assumed that one neutron is emitted per 
decay, the sum over all fission fragments of the decay 
rates given by Equation 4.17 is the rate of creation 


of neutrons at the model, 
~AS5 4) poe 
dn = >, ue (| e ) En, De _ (4.18) 


The absorption mean free path in the test section 
is about fifty times larger than the test section radius, 
therefore all neutrons freed in the test section will hit 
the model or walls. If all neutrons hitting the walls 


f= avsoOrbed, the ratio of neutrons hitting model to these 
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hitting the wall is given by the ratio of model radius to 
test section radius, Rm J Pe - For negligible 


model blockage of the test section assume this ratio is 0.l 


Then the delayed neutrons strike a unit area of 


the model at a rate 


277 Re Re art . (2s 19) 


The fission products yield are given in Reference 22. 


= 
ad 


An examination of the distribution of fission products 
shows that for a mass number less than 80 or greater than 
160, Ya is less than one part in tén thousand and the 
@Beneribution to Da is negligible regardless of the half 
life. Also, if ty, is greater than one minute, Aa is less 
meian .Ol1 See this combined with a maximum Ya of 
ten percent from Reference 22, shows that the contribution 
of any element with a half life much greater than one minute 
may be ignored in approximating the neutron flux due to 
secondary emission from fission fragments. References20 and 
22 were examined for fission products with yield fractions 
greater than on and half lives less than one minute. The 
isotopes satisfying these criteria are given in Table l. With 
these values, the rate of creation of neutrons at the model is, 
by Equation 4.18, 1.72 x 10°? neutrons per second and the 
neutron flux at the model due to secondary emission of neutrons 
is 


4 
Das = 2.3%*\0 De (4.20) 


It is noted that this method gives a conservatively 


high estimate of flux at the model due to assuming that 
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the total fission yield for a given atomic number indicated 
in Reference Bis the element with the atomic number with a 
short half-life, and further that these radionuclides all 


emit a thermal neutron upon decay. 


There are two unconservative assumptions in this 
development, that is that the initial condition is the 
absence of fission fragments at the beginning of the reactor 
Gwell time, and that only short lived neutron emitters are 
involved. In operating as a closed cycle,the supply of 
nitrogen and uranium will be contaminated by the fission 
fragments generated during previous residence in the reactor 
core, increasing the initial concentration of fragments and 


allowing build up of fission fragments with a long half-life. 


Due to the shear magnitude of the power involved in 
this facility, it is expected that the time for a given mass 
of fluid to complete a cycle will be long enough to allow 
decay in the concentration of the prompt neutron emitters 
and Beat provision can be made for monitoring and removal 
of long lifed fission fragments as a part of the fuel handling 


process. 


With these assumptions, the flux due to secondary 
emission from fission fragments is several orders of magnitude 
smaller than the direct streaming flux, and the former is 


neglected. 


The flux striking the model is then approximated by 


the direct streaming neutron flux of Equation 4.4 . 
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The average flux in the gaseous core is set by the 
reactor size, the power required, the energy per fission, é€, 
and the macroscopic cross section for fission of the fuel in 


the reactor core, 


PAT. | 
D ae E | | 


At a power level of 7100 megawatts with a 1.7 m core, 


the average flux in the gas core is 3.26 x Ome neuecen ene sec. 


Equation 4.4 gives the flux at the model-for these 


Gonaitions, 2.93 x 105° neutrons/ ane eee 


-~ 


Zeb Model Contamination 


The radiation hazard from the model is assumed to consist 
of secondary emission from neutron flux induced radionuclides 
on the model. The governing equation is of the same form 
as the build up of radionuclides in the core, then decaying 
after exposure. The build up of species A on the model during 


a test run of time,  . ,; as 


A(t,) = Za ee (1- e") (4.22) 


where Bs yp is the neutron macroscopic absorption cross section 
in the reaction ®Bten'—»A and B is one of the original 
materials of the model or a beans mute species. The rate of 


decay of A after the test is 


dA “rate, “4% 
oh - Au Bn (I-8*") (4.23) 


where t is the time since exposure to the neutron. After 


the decay rate for a species has been determined, the type 
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of decay must be examined to determine the mode of decay, that 
is alpha, beta, gamma Or neutron emission and the energy of the 
_ emitted particle. The particle type and energy set the maximum 
decay rate that a man can be exposed to and not exceed the 
limits of Ce rem/day. It is noted chee 1£ more than one 
species is formed or more than one mode of decay is present, 

the sum of the exposures will determine the allowable decay 


rates of each species or mode. 


As graphite is the prime candidate for material to be 


tested, it will be examined. Carbon undergoes the reactions 


Piel + en = oc Stable - (4.24a) 
ee a oc" =e tug. =5570 years (4.24b) 
ae + on = 4 + Ne stable (224) 
Or Ne. = 8.55 x mee ALOR Ene 

cae = 3 oe Toman ona 

Gers cia ae 9X oa oa 

te = 200 seconds 

Bove. = (ih ek (4.25) 


The initial decay rate for the model is, by Equation 4.23 


Aw _ gaxnio” GISMevh 
at ¢m® Sec 


With no absorption in the model, the surface beta flux 


for a cylinderical model is conservatively given by 


d Cus Vou = ol Cia Ry 
dt Au ou t e (4.26) 


Ror a carbon model with a radius of 1.25 centimeters 


df, - 5.27x lo Cis mex )/3 


dt Cm* sec 


cy 
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Reference 20 gives the exposure of betas of various 
energies which produces an allowable dose of .006 roentegen 
equivalent man per eight hour day. The exposure rate 
above is on the safe side, therefore a graphite model may 
be handled immediately after a run at Loo Megawatts 


for 200 seconds. 


Correlations between nuclear decay products and allowable 
dose rates are readily available in the literature an : 
A very slight protective layer, such as paper or CLOEh aivanel: 
provide protection against alpha particles, which do Woe 


Sa) 
Bimesent an external radiation hazard. 7 7 


It 1s also probable that high percentages of hydrogen 2 
will be present in an ablative material, and the decay of 
hydrogen isotopes yields a low energy beta particle in 
concentrations, determined by Equation 4.22, five orders of 
Magnitude below carbon and therefore,as a graphite model, 


also may be handled immediately after a test. 


This analysis can not possibly determine if each and 
every Material will be safe to handle after a test. The 
user will have to examine the reaction of each model 
individually to determine the contamination from each material 
and the time history of resulting total model contamination, 
as given by Raaeson 422 summed over all model components. 
The resulting model activity can be compared to allowable 


dose rates to determine if and when the model can be handled. 


ft is noted by examination of Equations 4.4, and 4.22 
fiat the decay rate at the model is directly proportional to 


the flux rate in the reactor. This in turn, by Equation 4.20, 
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ls dependent on the specific power of the reactor and varies 
inversely with the concentration of uranium in the core, so 
that low power density and high uranium concentration are 


desired. 


As is shown in Section 6.1, the operating power of 
the reactor is proportional to the area of the core, with 
core temperature and pressure fixed. The specific power varies 
then inversely as the core radius, as does the concentration 


of uranium in the core at small radii, shown in Figure 7. 


It can be concluded that, at least, for many possible 
re-entry materials there would be no difficulty in handling 
the models after testing and that model contamination is not 
a critical parameter in the design of this system, but should 
a particular material to be tested become highly contaminated, 
the degree of contamination can be reduced by reducing either 
the test duration or the specific power density. It is assumed 


that this does not occur in this study. 
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CHAPTER V 
REACTOR COOLING 
5.1 Heat Loads 


The walls of the’ gaseous Core reactor are subjected to 
Magh heat loads due to convective, conductive and radiation 
heat transfer and heat producing nuclear reactions in the 


Solid materials of the core materials. 


At the high temperatures of the core the heat load due 
to radiation transfer is much greater than the convective 
and conductive heat loads: the convective and conductive 
heat loads are considered negligible in this analysis. The 
core is’ assumed to radiate as a black body at 10,000°K and 


the walls absorptivity set as 0.6. 


The heat load due to nuclear reactions is assumed to 


be due to neutron and gamma fluxes from the core and is set 


; o ee 
eight percent of the power output divided by the core area. 


The heat load per unit wall area is 


Qe. = € Tle + .08 P/ Ae (5.1) 


where 
€ is the wall absorptivity, and 


@ is the Stephan-Boltzmann constant 


as 








42 
For the conditions of this study, the heaumleadsto the 


wall, given by Equation 5-1 is 


Qo, = 8.48 x 1o” col/o 2 Sé€c. 


pee Criticality Effects 


The loads given by Equation 5.1 are much too high for 
the thick beryllium moderator to conduct away from the core, 
and are predominantly due to radiation. Reference 7 suggests 
the use of a thermal barrier, cooled by transpiration 

and radiation seeding to protect the walls of a gaseous core 
nuclear rocket under similar high radiation heat loads. 

This method of cooling is not appropriate for use in the 
reactor used in this system. The large size of the reactor 
1s due to the high absorption cross section of the working 
gas,nitrogen. This causes a low value of the infinite 
multiplication factor in the reactor, Keo , given by 
Bguation 3.9. Examination of Equation 3-9 shows that 
increasing the macroscopic absorption cross section in the 
core will reduce the infinite multiplication factor, Kio. 

As this is already very close to 1.0 in this reactor, addition 
of even a small thermal barrier causes the reactor to be 


subcritical for any reactor size. 


Even if a material, with sufficiently low neutron 
absorption characteristics and high strength could be found, 
the transpiration flow necessary to internally cool the 
thermal barrier is not available due to the large area of 


this core and the relatively small throat area. 
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The moderator in this study 1s proceecreameneneene 
high thermal radiation heat loads by a seeded layer of carbon 
particles in the core near the walls. The small amount of 
“carbon necessary to protect the walls does not preceptibly 
lower the multiplication factor, but has a more pronounced 
effect on lowering the diffusion length, defined by Equation 
5.8 , sO that the carbon acts as a moderator and slightly 
decreases the critical size, as shown by Equation 3-10 ., 
This effect is due to the very small absorption cross section 
of carbon, which effects Kio , and the larger scattering | 
cross section, which reduces L . The effect of the carbon 


seed on criticality is ignored in this analysis. 
5.3 Reactor Flow Pattern 


In order to protect the moderator of the core, it is 
necessary to establish the following flow pattern. Consider 
any diagonal plane which includes the axis of symmetry of 
the nozzle. The reactor core flow and the nozzle are axi- 
symmetric. The reference direction,9=O, is aligned with 
the nozzle axis; the flow has a source at 6=fi , the fuel 
and nitrogen inlet, and a sink at the nozzle. A boundary 
layer is formed along the walls as the flow proceeds from 
inlet to the nozzle, much like a low speed divergent-con- 


vergent nozzle. 


The uranium fuel,in the form of a fine powder mixed 
with nitrogen, is sprayed into the core at the 9=T7 
position. The nitrogen is sprayed into the core through a 
annulus around the ewe nozzle. Perfect mixing of the fuel 
and nitrogen is assumed. It is assumed that the carbon seed 


tends to remain toward the walls and is limited to a layer 
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one-tenth of the core radius thick and starting at the core 


wall. 


The gamma and neutron fluxes to the moderator is un- 
aeftected by the thermal radiation barniaei eee eee a 
load to the moderator is the sum of the nuclear heat load 
and the thermal radiation which is not absorbed by the seeded 
layer. The moderator is cooled by a combination of radial 
injection of a fraction of nitrogen flow into the core, and 
water cooling. The nitrogen flow itself is not sufficient 
to cool the moderator and must also be used to establish the 


desired flow pattern. 
5.4 Radiation Barrier 


With the flow pattern in the core as established in 
paragraph 5.3, the flow is essentially flow in a nozzle 
with a diverging section proceeding the converging-diverging 
nozzle. Reference 12 gives the results of using a seeded 
layer to protect a gaseous core rocket nozzle. With perfect 
mixing of the fuel and working fluid, a seeded layer 1/10 
of the local nozzle radius thick will reduce the radiant 
energy flux by two orders of magnitude if the seeded layer 


absorption coefficient, K , is 1/3 ee eee 


Reference 13 gives a correlation ene meoese absorption 
coefficient per seed particle concentration, K/™, and the 
particle radius. Using carbon particles as the seeding 
material, it is desired that the average concentration of 
carbon wie les be minimized, both for nuclear and thermo- 


dynamic reasons. 
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The average carbon density, p , in the seeded layer is 


ai oe \ Ses = s a 


where 


Ps is the mass density of the carbon seed 


and Rs is the average seed radius 


. . =! 
With K set at VY, Cm,the average density is minimized 
ae ; 
at 3.7 £10 Jram/centimeter™ - this. occurs witha concenenactem 
10 ; 
of 435x100 eee ere having an average diameter of 


ol) microns. 


As the seeded layer does not fill the core, the average 


seed density over the core is 


6 An (RE -.9R2) (5.3) 
si Re 


f. - 


i 


a x O° 7 rams Jcentimete S 


The contribution to the macroscopic absorption cross 


section due to the carbon atoms in the core is given by 


Go 


teas 
£4: M,. 


Bifere cA is Avogadros number ancei ames the molecular weight 


= «5.4) 


Of carbon- 


For the concentration of carbon in the core due to the 
seeded thermal barrier, Equation 5.4yields an absorption 
: ° wt 
cross section contribution of \.6x*\0 a centimeter . AS 


this is four orders of magnitude less than the absorption 
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cross section of the nitrogen and uranium mixture in the 
@eore, the presence of the -carbon seed ,enmemirtca lit aus 


ignored. 


peo Moderator Cooling 


The seeded radiation barrier discussed in Section 5.4 
reduces the radiant heat load to the wall, given in Section 
Bil, by two orders of magnitude, but does not effect the nu- 
elear. radiation induced heat load. The resulting heat load 


to the moderator is 350 eave sec... 


The thickness of the moderator and the thermal conductivity 
of beryllium prevent this heat load from being conducted 


through the moderator and removed by back side cooling. 


A portion of the incoming nitrogen flow is radially 
passed through the moderator to absorb the heat load Re) tine 
moderator. The coolant nitrogen flow is injected tangentially 
into the core in such a manner as to be parallel to the flow 
of the boundary layer at the point of injection. This is 
to prevent the injected flow from driving the seeded thermal 
barrier away from the wall. A small amount of carbon, as 
is necessary to replace carbon in the thermal barrier which 
vaporizes or diffuses out of the thermal barrier, is carried 


by the cooling nitrogen flow. 


The maximum temperature of the beryllium moderator is 
set at 1000°K by structural strength limits-With one half 
of the primary nitrogen flow used to establish the desired 
flow pattern, the remaining primary nitrogen flow is capable 


of absorbing one third of heat load in the moderator, increasing 





AT 


in temperature from 293°K i 1000°K. The remaining heat 
load is absorbed by bleeding four percent of the water flow 
enroute to the water diffusor and allowing its temperature 


: O i 
to increase to 1000 K in the moderator. 
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CHAPTER VI 


NOZZLE HEAT TRANSFER 
6.1 Nozzle Size 


The hot gases in the core are radiating energy away 
from the core. As the moderator walls are not capable of 
transferring this heat load, they must be protected by a 
seeded thermal barrier, and the absorbed energy carried out 
of the core through the nozzle. Thus the minimum operating 
power is set by the radiant heat load from ENeS@EOTe, »ancesene 
Size of the core wall surface. The nozzle throat area is 
found by equating the radiative power in the corer comene 
energy flow through the nozzle, then, 


_ Qe. Ac 
A, = arm (6.1) 


For the operating conditions set in Section 2.2, the 


minimum throat size is 93 eantumecerae. 
6.2 Heat Loads 


Large heat loads are experienced by the walls of the 
nozzle as the flow is accelerated from stagnation conditions 
of the core to the high Mach number desired in the test 
section. The heat load is due to convective and radiative 
heat transfer from the hot gases and nuclear heating of the 


nozZle walls by alpha, beta and gamma interaction. To 
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prevent fallure of the nozzle, this heat leadeauctteepe 


transferred through, absorbed by or blocked from the nozzle 


Walls. 


Turbulent flow is assumed in the nozzle, also it is 
assumed that the flow is frozen,isentropic and one dimensional 
in determining local flow properties. The heat generated in 
the nozzle by nuclear interaction is assumed negligible compared 


to the convective and radiation heat loads. 


Sees Nozzle Configuration 


The nozzle consists of a series of longitudinal porous 
tubes, the area ratio of which is adjusted to provide an 
over pressure, of the cooling fluid in the nozzle tubes 
above the pressure of the flow in the nozzle, of 10 atmospheres. 
This provides the driving potential for transpiration cooling 
of the nozzle. Flow inside the coolant tubes is assumed to 


be turbulent. 


The tube cross section is a "U" shape, the flat sides 
of which may vary to provide the desired area ratio and also to 
facilitate joining the tubes together to form a smooth nozzle 
shape. The heat transfer surface is approximated by a circle 
whose radius is R. The fin effect of the radial sides of 
the cooling tubes is not accounted for, giving a conservative 


value for the heat transfer to the coolant. 


The coolant tubes support the thermal stress due to 
heat transfer across their radially inner-most face. They 
also support the hoop stress due to the coolant pressure. 


The tubes are closed by a shell which supports the pressure 
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of the hot gas in the nozzle. 
6.4 Maximum Heat Transfer 


_ The maximum heat load that can be transferred through 
the wall is set by the material properties, the configuration 


of the wall and the pressure the wall must support. 


For a thin pales tube, with’ the Wall thickness,d, very 
jweh less than the tube radius, Ry, the temperature gradient 
across the tube is linear. The thermal stress across the 
tube, in the tangential direction, are maximum in tension at 


the cool inner fiber, 


= +ae\ (T-. f2 (6.2a) 


and maximum in compression at the hot outer fiber 


Ow 2 -AE| (Te ez (6.2b) 


where 
M is the coefficient of linear expansion, 


E is the modulus of elasticity, 
lw is the hot side wall temperature, 


Tu.1ls the cold side wall temperature, 


‘and both & and E are functions of temperature. 


-The pressure differential across the cooling tube wall, 


p, causes a uniform hoop stress in the wall, 


s = + PR 
cl (6.3) 


where 


dis the wall thickness, 
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and the positive sign indicates tensile hoop stress, the 


negative compressive. 


The total stress is given by the sum of the thermal 


and hoop stress; scaler sum in this case as both are tangential. 


At the extreme fibers, 


GC = eee) a eS . . (6.4) 


The configuration chosen in this analysis utilizes a series 

of small cooling tubes forming the nozzle. This gives a 
tensile hoop stress and the maximum tobe at eaeee tensile at 
the cool inner fiber. The hoop stress relieves the compressive 
thermal stress at the hotter inner wall, where the material is 


weakest. 


The heat transferred through the wall is given by, 


a | 
@) as d (Ts - Te ] (6.5) 
where 


k is the wall thermal conductivity 


Equation 6.4 is solved for the temperature across the 


wall, and this is substituted into Equation G2. Su, 
, - PRe 
a 


KE a 
2k 


The maximum heat transfer through the wall, with respect to 


(6.6) 


Q= 


wall thickness,d, is obtained by setting, AQ Aid equal zero 





Ak? Ry _ 20k _ Oo 


aE d> % E <i? “(en 





— 
— 
a a. 


oe 


The condition for maximum heat transfer through the 


wall is, 


= Gay (6.3) 


z| 


G& _ OR _ 
Za a 


The maximum heat flow through the wall is 


- Sk . 
Cees a ZK EpRe (6.9) 


Here failure tccurs when the tensile stress at the inner 
fiber of the coolant tubes reaches the ultimate stress of the 
tube material. Radial stress is due to the pressure differen- 
tial across the tube See eee is negligible compared to the 
tangential stress. The nozzle has incroporated an expansion- 
slip joint at the up stream, low Mach end; this prevents 
longitudinal thermal stress so that the stress in the coolant 


mies is considered uniaxial. 


6.5 Material Selection 


To k 
A material was sought with a maximum value of a 





in order to maximize heat transfer with respect to material 
selection. This parameter is a function of temperature in a 
given material and the temperature must be high enough so 
that super-cooling is not necessary at the back side to absorb 
the heat load. Wall temperature also affects the heat load 

to the wall but examination of the Bartz's Equation on 
utilized to determine the convective heat load, shows only 


a weak dependence of convective heat transfer on wall temperature. 


a ; 4, oO 
This investigation chooses Rene 41° at give OL LOZa aK 


as having the desired structural properties. 
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At this temperature, 
Goris 146-7000 pounds Aen 
E is 2.4 x 10/ ound eee 
K is 0.05 cal/oem-sec- K 
% 1.57 x ind cm/em>K 


Minimum wall thickness is set by manufacturing limitations 


and is assumed to be 0.01 cm for this analysis. 


This choice of material and wall thickness gives a 


maximum heat transfer per unit wall area of 


| 3 
Quax = 194 410° dat /em’ sec 


provided that the cooling tube radius, set by Equation 
6.8 with a ten atmosphere pressure differential across the 
wall, is 5 cm or less. The temperature drop through the 


wall is 

Tn - Twe = 389°k 
— a = 636°k 
6.6 Test Section Mach Number 


The nozzle shape considered is a convergent-divergent 
nozzle formed by ae previously defined cooling tubes. | 
The nozzle is axa IM symmetric with a 45 degree convergent 
section, 15 degree divergent section and a radius of curvature 


at the throat equal to the throat inner radius. 


The heat load to the nozzle is weakly dependent upon 
throat radius, affecting the convecting heat transfer co- 


efficient only in this analysis. The pump work required of 
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the system as a whole, however, varies directly with the 


mass flow, which varies as the) area oF the tmucoane 


The throat radius is chosen in Section 6.1 as 5.45 
centimeters. 

The test section Mach number is arbitrary with a friction- 
less isentropic nozzle. In this nozzle, transpiration cooling 
1s required until the flow has expanded to about Mach 6, 

Pmereh requires a lengthy nozzle, of the order atm umeter—. 
Expanding the flow to Mach 6 increases the transpiration 
coolant requirement and nozzle losses. For this analysis, a 
test section Mach number of 3.0 is chosen to enable simulation 
of the desired flow fields, but reduce transpiration pumping 


requirements and nozzle losses. 


6.7 Transpiration Cooling Requirements 


The heat load per unit area to the wall, neglecting 


gamma heating, is given by, 


Q = QetQe (6.10) 


where Qo. is the convective heat load in the absence of 
transpiration .coolime, 


Qe is the radiative heat load, assumed unaffected 


by transpiration cooling. 


10 
The convective heat load was obtained by the method of Bartz 


for turbulent compressible flow, 


Qe =h( Tan -Ta) (6.11a) 


lB (BS). (BST EY TS sa 
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le (ie Steg]? [eset] F ere 


is 
is 
is 
is 
is 


is 


is 
is 
is 
is 
is 


1S 


the heat transfer coefficient, 
OpsG: 
the throat diameter, inches, 
the viscosity, 
the specific heat at constant pressure, 
the temperature, °R : 
sub = wall, 
Ww ‘ 
sub = bulk, 
Oo 
sub = adiabatic wall, 
aw 


the chamber pressure, 

the radius of curvature of throat, inches, 
l for the nozzle, 

the nozzle area ratio, 

the ratio of specific heats, 


the Mach number. 


: , 14 
The thermodynamic properties of nitrogen were used in 


this analysis, and the presence of the uranium was ignored. 


The radiation heat load is obtained by considering the 


gas to radiate at stream ambient temperature as a black body 


and the wall to have an absorptivity of 0.6, 


Q), = é. CT, ng KY en oT. (6.12) 


The radiation from the wall is negligible compared to the 


load to the wall and this reduces to 


Gar af in (6.13) 
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where 
€., is the wall absorptivity, 
€,, is the gas absorptivity, 
G is the Stefan-Boltzmann constant, 
T is the Temperature, ambient free stream “k, 


Woy is ene wall emissivity. 


The heat load to the nozzle wall in the absence of 
transpiration cooling is given in Figure 8 - The maximum 
heat load that the wall is capable of withseeming is lower 
than the load experienced. Transpiration cooling by the 
injection of low temperature nitrogen is used to reduce the 


convective heat load to the wall. 


The effects of transpiration cooling in reducing the 
convective component of heat transfer is given in Referencell. 
The Pecults are given in terms of the Stanton number reduction 
obtained by injection of a homogeneous coolant into a turbulent 
nozzle flow verses a blowing parameter, F/Ste. Here F is de- 
fined as the ratio of the transpirant mass velocity to the 


local free stream mass velocity, 
(2 G 
= ev (6.14) 
where G is defined as the transpiration mass velocity. The 
Stanton number is defined as, 
(6.15) 
eVicp 





S; = 


Then the Stanton number Gecter ied 1575 COscEnemrlrSte. Order, 


the reduction in convective heat load, 


S+/Sro = Ge / Gee (6.16) 
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The maximum heat load transferred through the wall, 
with transpiration cooling, is the sum of the reduced convective 
heat load and the radiation heat load minus the heat absorbed 


by the coolant passing through the wall. 


Qmee = Qu + Qy- Ge, (-T.) (6.17) 


where Cp is evaluated at local pressure and the average 


temperature in the wall. 


The convective heat load without transpiration cooling 


1s given by, 


Qe. = bh. (ie. -T } : (6.18) 
and _ 
h. Sto Cp " 


With Equations 6.18 and 6.19 the maximum heat flow through 


the wall can be rearranged to 


Qo — Qa. Ge _F C (h-7.) 
Q) oc Qe Qo. S;, Co ( Taunt) (6.20) 


The left hand side of enn equation is known at each 
station and the right hand side suggests the possibility of 
a graphical solution using the difference between the Stanton 
number reduction data of Reference 11 and a family of Seeeane 
lines of slope Cy (TW. -Twe) /Cp (Tau - TL, ) . Figure 9 
reproduces the Stanton number reduction data of Reference 11 
and illustrates the graphic method of solution. The blowing 
parameter, F, is obtained by this method and from it, the 


transpiration mass velocity necessary to reduce the heat 
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load through the wall to the maximum allowable. This flow 
is graphically integrated to obtain the total required 


transpiration flow. 


The transpiration cooling flow per unit axial nozzle 
length, 21TRG , 1s given in Figure 10 , along with the total 


transpiration mass flow. 


Under the conditions of this study, the transpiration 
coolant flow necessary to prevent thermal damage to the 
nozzle and test section is 1\1.25~x 10° qram / Second , and 
the ratio of transpiration flow to the primary flow in the 


nozzle is, 


ee ey 
wv} 


This estimate on transpiration cooling requirement is 
conservatively high as the flow is assumed to radiate as a 
black body and the effects of thermal radiation blockage by the 


seeded layer in the core is neglected. 


6.8 Secondary Nitrogen Flow ° 


The heat flow through the wall must be absorbed by 
the flow in the coolant tubes. The flow properties in the 
tubes are set by matching the heat transfer at the throat 


as this is the point of maximum heating load. 


The flow in the tubes can not be set point by point 
to enable maximum heat transfer to be possible along the 
entire length of the nozzle and heat transfer rates below 
way are possible along the nozzle. This problem is 
recognized but not investigated further in this analysis 


as it is felt that the transpiration cooling flow can be 
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me : ; . . 
tailored to fit a wall with variations in heat transfer 


capabilities. 


The heat flow through the walls is absorbed by the 


coolant flow, 


One = he (Te=Te) (6.21) 


Here radiation and nuclear heating of the coolant are ignored. 
-The coolant is flowing in a nozzle shaped duct but here 

fluid is bled out of the tube. In the absence of data on 
sucking of a turbulent boundary layer, the data of Reference 11 
is extrapolated to include the sucking case. This is also 
illustrated in Figure 10, and is analogous to the treatment 


of a sucked laminar boundary layer given in Reference 25. 


The method of Bartz is used as before to determine the 
heat transfer coefficient. The coolant tube nozzle throat 


meas is set at .5 cm. 


The chamber temperature and pressure of the secondary 
nitrogen flow is set by wetaelaling the heat transfer through 
the wall at the throat to the heat apaeebedune the coolant 
of that point, taking into account the effect of the transira- 
tion cooling flow sucking the boundary layer in the coolant 


tubes. 


For the conditions of this study, the secondary nitrogen 
chamber pressure is 1010 atmospheres and hs temperature is 
530K. With these chamber conditions the mass flow of 
secondary nitrogen is 2.17 x 10° grams per second and the 


ratio of secondary to primary nitrogen flow is 0.67 to l. 
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The secondary flow of nitrogen that is used to provide 
pee piteacnan and backside cooling of the nozzle is approx- 
imately the same as the primary flow through the nozzle. This 
doubles the amount of sane work required, bit as) che snozzile 
would fail without the transpiration protection, this additional 
pump work is an acceptable penalty that must be incorporated 


into the design of the system. 
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CHAPTER VII 


AUXILIARY EQUIPMENT 
7.1 General 


Operating this facility as a closed cycle removes 
some of the technological problems involved in the design 
of a gaseous core rocket. The problem of containing the 
gaseous fuel in the reactor core is avoided here by operating 
as a closed cycle. This introduces the problem of additional 
equipment in the cycle. Among the additional equipment needed 
are high powered turbines, high capacity pumps, large heat 


exchangers and fuel separation and processing equipment. 
7.2 Water Diffusor 


The flow at the exit of the test section is mixed 
with the secondary nitrogen flow which was utilized to provide 
Beenspiration and backside cooling of the nozzle. At this 
meme the flow stagnation temperature, determined from the 
enthalpy of the mixture, is 8000°K and the stagnation pressure, 
assuming isentropic expansion through the nozzle, is 1000 
atmospheres.Conventional heat exchangers can not operate 
at these temperatures and pressures, neither can conventional 
turbo-machinery. To lower the stream's temperature and pressure, . 
a water spray diffuser 1s employed. The water is sprayed 
directly into the hot nitrogen stream and is heated to super — 


heated steam by the nitrogen, which loses energy in the process. 
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The cross sectional area of the diffuser can be designed 


to obtain the desired outlet conditions. 


Conservation of energy determines the required water 


flow. Conservation of energy requires that 


mo Hu (Ta, R.,) ema H,. (To,R.) = wm He (Tor $B.) & mat. Mo, tI-$12.) 


where Y\ is the mass flow rate vo 
H is the total enthalpy 
5 is the mole fraction of water 
I-§ is the mole fraction of nitrogen 


and subscript w refers to water 


1 inlet conditions 
OQ stagnation conditions 
2 Ouclet conditione. 


Equation 7.1 was utilized to obtain the required mass 
flow of water necessary to reduce the temperature and pressure 


of the nitrogen stream to the desired values. 


The size of the water diffuser is set by the perfect 


gas law, written in the form 


> 


Ro. Am ™ 
PR, Ar ™, Mi 





‘Va ¥, Me 
The My 





(7.2) 


Where all values are static values and 


KH is the mixture average molecular weight 


ysis the mixture average specific heat ratio 


With the properties of nitrogen given in Reference 
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14 and those of water from Reference 26, Equations 7.1 
and 7.2 can be used to find the necessary water flow and 
@LEEUSOY Gxit sarea, Once (ENS vee eMache nun De yamec noord eins 


and pressure are specified. 


The diffusor outlet total temperature is set at 1000°x, 
aS an upper operating temperature for turbo-machinery. The 
outlet total pressure is set a 100 atmospheres and the 


Mach number at O.OL 


With these conditions set, the required mass flow of 
, 5 
water is 9 x 10 grams per second, and the output area 


Greene diffusor is 8.25 x one a 


fee Fuel Separator 


The cooling process in the water diffuser cools uranium 
Fuel in the flow well below its melting point. The super- 
cooled uranium rapidly condenses and forms small droplets 
and particles which are suspended in the flow. The thermal 
energy of the flow is used to provide a centrifugal force 
field to separate the particles from the gas. The drag on 
these particles is given for Stokes flow as 
& . 454V (7.3) 
m Re Pu 
where 4 is the viscosity of the steam, nitrogen mixture. 
The equation of radial motion of a uranium particle in 
@ecemtrifugal force field ae with the multiplier of velocity 


in Equation 7.3defined aoe 


r +Cor -wr =o (a) 


foes has a sOlution of the form 
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Vo _ Cor lor rAur , _ Can flor + 4k 
a 


= (7.5) 
C3. +4? C 


For small particles, the drag coefficient Cs is 
large and the second exponential term of Equation 7.5 
decays rapidly regardless of the value of @) and the 
first term determines the radial position of the 
particle. For small values of as , the particle is strongly 
restrained by the drag forces of the fluid and a centrifugal 
separator is not very effective. For a centifugal separator 
to efficiently separate the small particles, the rotational 
velocity must be at least of the order of the drag cofficient. 
Under the conditions of the flow just after the water diffusor, 
the angular velocity must be of the order of 100 radians 
per second for micron sized particles. This imposes large 
Velocity requirements if the flow is swirled in one section 
at the outlet of the water diffusor, so the flow is divided 


into several centrifugal separators. 


7.4 Turbines and Pumps 


The ideal power available in the flow as it comes out 


of the uranium separators is given by 
P= Mw H.. rm He (7.6) 


It is convenient however to expand the flow only until 

the water vapor in the flow reaches its saturation point. 
At Bhi s point, the water can be condensed and removed from 
the flow stream and the nitrogen further cooled to reduce 
the pump work required to increase its pressure to core 


entry conditions. Under these conditions Equation 7.6 
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becomes 


P= mi {Ha (rock feisoatn) — WulTeor Fe hae)] 


+ My {H (loos k, Licooakm ) = Ho Gus P, .« } (77) 


For isentropic expansion of the water vapor, the ideal 


power available is 2,520 megawatts. 


After condensation, the water is pumped as a liquid 
to 120 atmospheres and injected into the water diffusor. 


This requires an ideal pump work, 
: . 
me - ae, Sa (7.8) 


Ideally, 10.9 megawatts are required for the water pump 


work. 


The nitrogen is cooled in a heat exchanger to 293°K 
and then pumped back to 1000 atmospheres. This must be 
done in at least two pumps with interpump cooling to maintain 
the temperature at or below a structurally safe limit. The 


ideal nitrogen pump work is given for each pump by | 


We. = Mma H(t, 8.)- Hah, R ) 7.3) 


where the temperatures and pressures are related by the 


isentropic compression relation 


_ Cos 
Naess (B:.) S : (7.10) 
aYr He 
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The total pump work for all nitrogen pumps is 868 


megawatts. 


The total heat rejected in the cycle is the sum 
of the heat of condensation of the water vapor, the 
cooling of the nitrogen prior to being pumped, and inter- 
pump cooling after pumping so that the nitrogen can be 
used to cool the reactor moderator. The total heat rejected 


mm the cycle is 7.1 x 10° megawatts. 


bi Walter as usedgas a heat sink for this rejected heat, 


the mass flow required is 


Qre; 
Mere, * cae | (7.11) 


O 
For a temperature rise of this coolant water of 100 K, 


Eueemass flow required is 1.1 x Te grams per second. 
7.5 Minimum Power Level 


The minimum operating power level of this test facility 
is set by the criticality requirements of the reactor, Chapter 
III, and heat transfer considerations in the reactor and nozzle, 
Chapter V and VI. This high operating power requires high 
mass flow rates and large capacity pumps and turbines, as 


gmdicated in Section 7.4. 


There is no lack of power available to drive this machinery 
and small units can operate in parallel to provide the desired 
Capacity. As the system is proposed, there is a requirement 
that the combined efficiency of the turbines and pumps be 
about 35 percent, but if all of the energy in the cycle is 


used, efficiencies as low as 12 percent are acceptable. 





or, 


Equations 5.1 and 6.1 can be combined to express the 
Operating power minimized'’with respect to the area of the 


throat, 


P-(ects 4.08 P/A.) Ac (712) 


The area of the core is minimized for any pressure 
and temperature by setting the ratio of nitrogen mass 


flow to uranium mass flow at 15.4, as was shown in Section 3.2. 


Equation 3.11 shows that, with the mass flow ratio 
fixed, the critical radius, R, Varies anvercelly with the 
number density of uranium, Ny, . For constant operating 
pressure, the number density of uranium varies inversely 
wmen the Operating temperature, |, . The critical radius 
is then proportional to the operating temperature and the 
core area is proportional to the square of the operating 
temperature. 

With fixed mass flow ratio and operating pressure, the 


minimum operating power is, to the first order, 
Ey Cy (7.13) 
mn | © ° 
where C is a proportionality constant. 
Equation 7.13 is valid only in the temperature range 
where the heat loads to the moderator are too large to be 
removed by conventional cooling and the energy absorbed in 


protecting the moderator walls must pass through the nozzle. 


A plot of the estimated minimum operating power as a 


ee ee ee we AP APP ane om 
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function of temperature is given in Figure ll. At operating 
temperatures below 7000°K, this facility 1s in the power 
category of large arc tunnels and MHD devices, and the 
relative merits of each type of facility would navemrtoube 
considered before deciding which system is most desirable. 
At the higher temperatures at which this facility is con- 
ceptionally most useful, the minimum operating power rises 
rapidly, as well the size of the physical plant. It is also 
noted that reducing the Operating tempeeraicuremmedmecc Che 
re-entry simulation capability of the test facility, as 


illustrated in Figure 2. 


The economics of developing such a facility must be 
examined to determine if the need for such a testing capability 


justifies the expense of so large an undertaking. 
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CHAP aR V It 


CONCLUSIONS 


Within the scope of the material examined in this 
report it appears to be within the capability of projected 
technology to design and build a high enthalpy test facility 


powered by a gaseous core reactor. 


This test facility can not exactly duplicate all of 
the conditions imposed upon a vehicle during an atmospheric 
re-entry, but does provide a means of testing re-entry 
materials and shapes under the high temperature and shear 
conditions which are encountered during a re-entry. Existing 
facilities cannot attain the high temperature, high eee ee 
long duration flows necessary to test the shear mode of 
failure. The gaseous core reactor provides the power required 


for these flows. 


Nuclear contamination of several proposed re-entry 
materials is slight enough so that the model may be examined 
immediately after a test run. This may not be true in 
general and each proposed test material must be examined 
by the methods of Chapter IV to determine the degree of 


contamination. 


The high neutron absorption cross section of nitrogen 
as a working fluid causes the reactor to be large and the 


mass flow ratio of uranium to nitrogen to be relatively 





a 
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high. The reactor size and operating power is minimized 
‘by setting the mass ratio of nitrogen to uranium in the 
core at 15.4 to 1. This adversely affects the quality 
of the mixture for simulating atmospheric conditions, 


but is slight enough to be considered acceptable. 


The minimum operating power of the test facility is 
set by the radiant power of the core which cannot be 
absorbed by the moderator and must pass through the nozzle. 
In the neighborhood of 10,000 are for a constant mass ratio 
the minimum. operating power varies as the temperature in 
the core to the sixth power and as the inverse of the core 


pressure. 


‘Operation of the facility at the minimum power level 
requires high mass flow rates and associated high capacity 
pumps and turbines. These pumps must be capable of high 
pressure and high pressure ratio operation. Efficiency 
of the combined turbine and pump can be as low as twelve 


percent for self sustaining operation of the system. 


The turbo-machinery required set a practical upper 
limit on the operating power of the system. Structural 


considerations impose an upper pressure limit on the system. 


This analysis has examined only the minimum powered 
point for the operation at the maximum pressure and temperature 
and has found that such a facility appears technically feasible. 
Operation of a facility at lower power levels is possible 
with a reduction in chamber temperature and associated loss 


of re-entry simulation capability. 
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It is recommended that the economics:involved in the 
construction, maintenance and operation of this system be 
examined to determine if such a large facility would be 


economically practical. 


72 


TABLE I 


FISSION PRODUCTS 


Radionuclide Half Life Yield Fraction _ Decay Constant 


A | ti, (sec) Via Aa (sec) 
ee 55.6 4973 1.2 x 10 - 
pr °7 4.51 5 x 10> 1.54 x 10> 
Ag < Dae? nome 29 ae 
In oe qe jon 9-6 x 100° 
= ae 13 | 10-7 5.4 one 
eb = 1.52 4% 101- eee laa 
a! 22 ce One Se aaearon< 
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